Phase transition between the cholesteric and twist grain 

boundary C phases 



I. Luk'yanchuk 1,2 '* 
1 L.D. Landau Institute for Theoretical Physics, Moscow, Russia. 
2 Departamento de Fisica, Universidade Federal de Minas Gerais, 
Caixa Postal 702, 30161-970, Belo Horizonte, Minas Gerais, Brazil 

(February 1, 2008) 



Abstract 



The upper critical temperature T C 2 for the phase transition between 
the Cholesteric phase (N*) and the Twist Grain Boundary C phase with 
the layer inclination tilted to the pitch axis (TGBct) in thermotropic liquid 
crystals is determined by the mean field Chen-Lubensky approach. We show 
that the N*-TGBct phase transition is split in two with the appearance of 
either the TGBa or the TGBiq phase in a narrow temperature interval be- 
low T C 2 . The latter phase is novel in being superposed from two degenerate 
TGBct phases with different (left and right) layers inclinations to the pitch 
axis. 
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I. INTRODUCTION 



A Twist Grain Boundary (TGB) state that appears as an intermediate state at the 
Cholesteric (N*) - Smectic (Sm) phase transition in chiral thermotropic liquid crystals 
was predicted theoretically by Renn and Lubensky [Q in 1988 and then, one year later, 
was independently observed experimentally 0|§. Since that time a wealth of properties of 
this new state were discovered in a number of experimental 0-13] and theoretical 0-17 
investigations. 

The results of these studies and the results of the present paper are summarized in the phase 
diagram of Fig. 1 where the parameters t, <Jj_ (whose meanings will be explained in Sect. 
II) are controlled by the following experimental conditions: temperature, concentration, 
pressure etc. The reason for such a variety of intermediate phases is that, the direct N*- 
Sm transition can not occur in a continuous way since the cholesteric twist of the director, 
n (r) = (0, sin (k x) , cos (kox)), is incompatible with the smectic layered structure. The last 
one in chiral liquid crystals is known to be of either SmA or SmC* type: in SmA the 
director n is parallel to the layers modulation vector q whereas in SmC* it is tilted with 
respect to q by a constant angle 8q and forms a conical precession along the normal to layers. 
Therefore, the transition occurs either by the first order untwisting of n (r) or via formation 
of intermediate TGB phases. The actual sequence of the intermediate phases depends on 
the final Sm state that occurs at low temperature. There is only one intermediate phase 
(TGB a) when the transition goes to SmA and a series of phases when the transition 



goes to SmC* fl5] 



The general structure of the TGB state is shown in Fig. 2a. The compromise between the 
cholesteric twist of n (r) and layered structure of SmA, C* is achieved by formation of a set 
of rotated smectic slabs (blocks), normal to those being follow the pitch n (r). The slabs are 
separated by grain boundaries consisting of a series of equally spaced screw dislocations that 
provide the junction of the layers in adjusting slabs. The slab width £&, dislocation spacing 
Id, director pitch P, and layer spacing d are related by the following topological constraint 

2nl b l d = Pd (1) 

Coupling of the director with the modulation vector q results in the unbending of n (r) 
close to the block center. When temperature decreases the director pitch and slab width 
diverge with further untwisting of n (r). Finally, they tend to infinity corresponding to the 
transition to the Sm state. 

The variety of TGB phases in Fig.l is provided by the different internal structure of TGB 
blocks that are shown in Fig.2b-e. Generally, TGB blocks are reminiscent of the final Sm 
state that occurs at low temperature. Therefore, the TGBa block is just the SmA slab, 
shown in Fig. 2b, confined by grain boundaries, the layers modulation vector q being parallel 
to the director n in the slab center [I] . Similarly, the TGBq slab is provided by smectic layers 
that are inclined to n by the angle ~ 6q. The inclination, however, can be done in other 
ways: when tilted layers are either parallel to the pith axis x as in Fig. 2b, or tilted to it as 
in Fig.2c. We call these phases TGB Cp and TGB C t- The blocks of the TGB Cp and TGB C t 
phases actually have the structure of slabs of differently oriented SmC, as distinguished from 
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its chiral analog SmC* by the absence of the director precession. Under certain conditions 
a transition SmC — > SmC* occurs inside the TGBc slabs ||15|| . The corresponding TGBc* 
phase that appears close to the bulk SmC* phase (see Fig. 1) will be not considered in this 
paper. 

Originally the TGBc p phase was assumed to be an intermediate TGB state at the N*-SmC* 
transition Jl5] (there it was called TGBc)- However, x-ray experiments |7j and theoretical 
estimations []I7| demonstrate that the TGBct phase is indeed more stable. In |T7| this 



phase was called as the Melted Grain Boundary (MGB) phase to stress that the smectic 
order parameter vanishes at grain boundaries because of the small distance between screw 
dislocations. We prefer, however, to use the TGBct notation to emphasize the geometrical 
structure of this phase. 

In this paper we revise the calculation of for the upper critical temperature T& (M4- 
Mq-Mi-Mc line in Fig. 1 for the N*-TGBc transition, taking into account the recent proof 
of the stability of the TGBct phase [|T|,|iT| that was not considered in [T3|. We confirm and 



expand the estimation of JT5| to the whole region of parameters. In addition we calculate 
the principal parameters of the TGBct phase: upper critical temperature, slab width Z& and 
the X-ray diffraction pattern that can be measured experimentally. 
Several features that modify the phase diagram calculated in |TB|,|T3J follow from our anal- 
ysis. The TGB a phase that was shown in [|T|,|T4|| to be stable when the transition is to the 
SmA phase, penetrates also into the SmC* region as a narrow stripe in between the N* 
and TGBct phases, and finishes in the tetracritical point Mq far inside this region. The 
N*-TGBct transition is split either by this TGB a stripe or by the narrow region of the new 
TGB 2q phase. The TGB 2q slab is superposed from two equivalent SmC populations with 
left and right layers inclined to the pitch shown in Fig. 2e. 

This new phase can be viewed as a kind of standing density wave quantized by the grain 
boundaries. In reality it can be observed in between the N*-TGBA-TGBct-TGB2 q tetra- 
critical point Mq and the point M\. Location of these points is calculated in this paper. 
The upper critical temperature T c i should have a kink at Mq. Under certain conditions the 
enhancement (oscillation) of T c2 between M Q and Mi can be observed. 

A remarkable feature of the TGB a phase was noted in |I[ to be an analogy with the 
Abrikosov vortex state in superconductors in a magnetic field. This completed an analogy 
between the superconducting transition in metals and the phase transition between Nematic 



(N) and SmA phases in liquid crystals, first pointed by de Gennes (T^JT^|. In the present 



article we show that the TGBct phase is the analog of the mixed state in superconduc- 
tors with a space-modulated order parameter (like Larkin-Ovchinnikov-Fulde-Ferrel phases 
PD|PT[ ]), providing that modulation is perpendicular to the magnetic field. Another inter- 
esting analogy we discuss is the similarity between the TGBA-TGBct transition and the 
transition between the symmetry differing phases in an "unconventional" superconductor 
UPts in the magnetic field p2| . 



II. BASIC EQUATIONS 
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A. The chiral Chen-Lubensky Model 



On a quantitative level, the appearance of TGB state is described by the Chen-Lubensky 



(CL) model | 23fl , which is known to be a quite general approach in explaining various phase 
transitions between cholesteric (nematic) phases and modulated smectic phases. In this 
model the cholesteric and smectic phases are described by two coupled order parameters: 
by the twisted director n{r) and the space modulated complex function ip(r), where the 
modulation of the smectic density is given by the real part of ip(r). The resulting energy 
consists of two parts: 

J~CL — J~il> + J~F (2) 

where the elastic Frank energy: 

Tf = -K\ (div n) 2 + -K 2 (n • curl n — k Q ) 2 (3) 



2 v 7 2 
+ -i^ 3 (n x curl n 
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provides the twisted texture of the director. In the cholesteric phase n (r) = 
(0, sin (fco^) , cos (fco^))- The chirality is provided by the parameter k$. When ko = 0, 
expression (||) reduces to the elastic energy of the nematic phase with n = const. 
The smectic state is described by the Ginzburg-Landau (GL) functional J 7 ^, which has a 
fmite-g instability for the order parameter ijj(r) provided by a gauge derivative D = V — iqon: 

F i , = a(T - T NA ) | V | 2 I ^ I 4 n i n j ( 4 ) 
+C ± (5 i:j - mnj)) (Dip) (Dj0)* + D (D 2 ^) (DV)* 

We take the quartic in the gradient term in the isotropic form D (D 2 ^) (D 2 ^)* which 
is slightly different from the original CL model where this term was written as: 
D± (Sij — rurij) (5m — rikrii) (DjDj^) (DfcI}?/>)*. This does not change the final results but 
simplifies the calculations. 

Now we give a brief review of the properties of the CL model. In the nonchiral case 
(when k = 0), a second order phase transition from nematic to smectic phases takes place 
||23|| . The type of the smectic phase depends on the sign of C± (terms with Cm and D are 
assumed to be positive). When C± > 0, a transition to the SmA phase with the order 
parameter ip(r) ~ ip e tqonr occurs at T = T^a- When C± < an additional transversal to n 
modulation occurs, resulting in a SmC order parameter ip(r) ~ ipoe iq}r with q = (gon, q c ). 
This modulation is stabilized by the quartic gradient term at wave vector qc = (—C±/2D) l l 2 
which gives the layers inclination angle as #o = arctan(gc , /o , o) = arctan(— C ^jlDq^) 1 !' 1 . The 
N-SmC transition occurs at T^c = Tna + Cj_/4D. 

In a chiral case |14| (when ko ^ 0), the situation is more complicated since the gauge 
derivatives and anisotropic coefficients in the gradient terms in Eq.(^) depend on the direc- 
tion of n which is not uniform in a space. This provides a coupling of the order parameters 
n(r) and ip(r) leading to an untwisting of n(r) when the smectic layers are formed. Such 
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a process occurs either directly, by the first order transition, or via the intermediate TGB 
state. In both cases the final smectic phases are again either SmA or SmC. 
The scenario of the phase transition depends on the parameters of the CL model. We 
formulate those conditions below. When the direct N*-SmA, C* transition takes place, the 
critical temperatures are calculated by comparison of the energies of the N* and SmA, C* 



phases | 14}j . The N*-SmA-SmC* tricritical point in the (T, C±) plane was found to be: 

(T*, C* ± ) = (T NA - (jgK 2 )V 2 ko/a, (gK^k^aK^) (5) 

The transition lines are given by: 

T N . A = T*, C ± > CI (6) 
Tag* — Tna + {Tn*a — T NA ) ■ C* / C±, C± < C]_ 
T N * C * = T N * A + ( C ± - Clf/AaD, C ± < C* ± 

When the transition occurs via an intermediate TGB state, the critical temperatures and 
the detailed structure of the TGB state are provided by a nonuniform solution of the GL 
equation obtained from (|4|): 

ip + g^U (7) 




nD)> - D \-jy + ^ j i/> 
where, A 2 ip = A (Ai/S^ . 

In this paper we are interested in the structure of the TGB state just below the N*-TGB 
phase transition that takes place at the upper critical temperature T&. We consider the case 
when C± < 0, that is when the system has a tendency to form the SmC* phase at low 
temperature. In the following subsection we derive the basic equations that describe this 
transition. 



B. Cholesteric-TGB transition 

The order parameter in the TGB state below T C 2 is given by a periodic superposition of 
rotating blocks equally shifted by a distance lb 

^ (r) = ]>V 7m /(x + ml h )e l ^ mr (8) 

m 

The component of the modulation vector 

transverse to the pitch q_|_ m = q (0, sin (komlb + <p) , cos (k mlb + y?)), follows the director 
twist: n(r) = (0, sin (k x) , cos (k x)) with phase advance (retardation) (p. Factors e* 7m are 
the degenerate phason degrees of freedom [|l|]. 

The block profile function f(x) is localized within the block width If, and, together with the 
phase ip, provides the structure of the TGB slab. In the TGB A phase (p = and f(x) is a 
centered bell-shape function (which has a Gaussian profile exp(— k x 2 /2q ) at C±_ 3> 



5 



The structure of the TGBc p slab (Fig. 2c) is given by a centered bell-shape function fix) 
and nonzero advancing (retarding) phase angle ip ~ 6 Q [15 that results in the tilting of the 
layers parallel to the pitch axis . In the TGBct phase the angle ip is equal to zero and f(x) 
is modulated in the x direction as exp(±i9x). This corresponds to the right (left) transverse 
inclination of the layers as shown in Fig. 2d. 

The profile function f(x) of the slab, located at the origin, and phase <p are found by solving 
of Eq.QTp with the substitution tp (r) = f(x)e ici±or ~ f(x)e iqo{z+vy) . 

Several simplifications are used in Renn-Lubensky theory. Usually the block width 1^ is 
much smaller than the cholesteric pitch P = l-njk^. Therefore, on the scale of Z& the twist 
of n (r) is minimal and approximately is written as n (r) pa (0, k x, 1 - (k x) 2 /2). Next, 
just below T c2 the amplitude of f(x) is small and only the terms linear in f(x) are relevant 
in Eq.(0). 

The corresponding linearized equation for / (x) and ip is: 



a(T- T NA - C 2 JAD) f = -Hf where (9) 



H^- D(-dl + {q k,)\x- V /k,f + ^ 



+ (C l \-C L )^(x 2 -2xp/k f 

This equation has a set of localized eigenstates f n (x) with a discrete spectrum of eigen- 
temperatures T n . The upper critical temperature T c2 of the N*-TGB transition is provided 
by the maximal value of T n . The block profile function is given by the corresponding eigen- 
function f nc2 {x). 



We assume further that <p = 0, that is, the parallel layers' inclination according to [|lj] does 
not occur. The justification of this approximation will be presented in Sect. HID. 
It is convenient to use the dimensionless units: 

x = q x, d x = q^d x , a\\ = — L a ± = —^j*, (10) 

t = a(T- T NA ) /Dql b = k /q , H = H/Dq* 

We discuss first the values of parameters b,a±,a«. Parameter b has a sense of the power 
of the cholesteric twist. Usually in chiral liquid crystals the interlayer smectic spacing is 
much smaller than cholesteric pitch having b 1. The small tilting angle 8 results in 
<j± = tan 2 $o < 1. In contrast, parameter cry, which is related to the ratio of the layer 
compression elastic constant Cnq^ to the layer curvature energy Dq$, is larger than one 



T7,2|] and therefore is much large than a^. 

Neglecting C± in comparison with C\\ in the last term of Eq.([|), we rewrite it in dimen- 
sionless units as: 



(t - a \)f = -Hf where (11) 
H = (-d 2 x + b 2 x 2 -a ± f + ayfeV 
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The order parameter tp (r) of the TGB state is reproduced by substitution of corresponding 
eigenfunction f(x) into (|]). Then, one can calculate the free energy Tql (1) of the TGB 
state as a function of the slab width Z& or, more conveniently, as a function of the geometrical 
factor hi Id- The actual value of this ratio is found by minimization of Tcl with respect to 
lb I Id- It was found [1[] that h/h — 0.9e when C± ^> and the TGB a phase has a Gaussian 
profile in the slab. The factor e depends the relative strengths of cholesteric splay (bend) and 
the twist elastic energies Kx^/K 2 - It varies from e = when K\^/K 2 = (the calculation 
in this case is reduced to minimizating the Abrikosov factor (3 = (ip A ) / (tp 2 ) ), to e ~ 1.5 
when K\^jK 2 is large. We will use this result to evaluate h/ld in the TGBA,ct phases when 
C ± < 0. " 

III. RESULTS 
A. General 

In this section we calculate the eigenstates of Eq. ( |TTD that correspond to the upper critical 
temperature t c2 , and discuss the block structure just below t c2 for the different types of TGB 
state at a± > 0. 

Note that eigentemperatures t n (including t c2 ) and eigenfunctions f n (x) are generally func- 
tions of b,aj_,au. Because of the scaling properties resulting from Eq.flTl|): 

£„(6,cr||,ojJ = o\ ■ t n (b/a±,a\\, 1), (12) 
f n (x,b,cr\\,a±) = f n (ay 2 x,b/a±,a\\, 1). 

parameter <t_i_ can be excluded if one considers the dependence t n (b) on the rescaled coordi- 
nates b/a±,t/crj_. 

To proceede with the diagonalization of Eq.flTTI) for <r_|_ > 0, recall first the results of P|JT^] 
for the opposite case of <r± < when a transition occurs to the TGB a phase. When |<t_l| is 
large enough, the operator 7i was shown [|TJ to be truncatable to the more simple form of 
the harmonic oscillator — 2a± (—d x + b 2 x 2 ^j + a\. The lowest eigenstate gives the Gaussian 

profile e~ bx2 of the TGB a block and the upper critical temperature t c2 = 2a±b. This result 
can be improved if one considers the residual part of H as a perturbation: 

t c2 = 2a±b- (1 + 3(T||/4)6 2 (13) 

The second term is smaller than the first one if —<J±/b > 0.5 + 0.38<7||, which is the condition 
of applicability of the approximation. The opposite case of a±/b ~ will be discussed later. 

Figure 3 shows the result of numerical diagonalization of (|TTD for different an when a± > 0. 
We use the t/aj_, b/a± coordinates to trace the two highest eigentemperatures (the maximal 
one corresponds to t c2 ) as function of b. Note first the concurrence of two eigenstates with 
close eigentemperatures t + and t_ that give the upper critical temperature t c2 = max(t + , £_). 
These levels often cross when b changes, resulting in the oscillating behavior of t c2 (b). The 
oscillations of t c2 {b) are strong at cry = and more pronounced still at large values of a\\. 

The appearance of these nearly degenerate eigenstates is related to two equivalent, right 
and left inclinations of layers in the TGBct slab. The corresponding block profile functions, 
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}'r{x) and however, can not be the eigenfunctions of Eq.flllD since they do not possess 

the definite parity with respect to the symmetry operation x — > —x, which is a property of 
H. The proper eigenstates f± with eigentemperatures t± are constructed as superpositions 
of both populations as f± = /t that results in the TGB 2q phase just below t c2 . The 
function f±(x) has cos- (sin- ) like oscillations and can be viewed as a standing wave between 
the grain boundaries. The effects of commensurability of this wave with block width lb 
stabilize either the even cos-like or the odd sin-like behavior as the lowest state to have a 
vanishing order parameter at the grain boundaries. The block width 4 changes as a function 
of b that alternates the order of the f± eigenstates leading to the oscillations in t C 2(b). 
The tendency to SmC slab formation at lower temperatures results in the further transition 
from the TGB2 q to TGBct phase determined by the lowest from t+, t_ critical temperatures 
renormalized by the nonlinear term g \ip\ 2 ip in Eq.(|7]). Below this second transition, one of 
the populations starts to be suppressed and a block profile function is constructed now from 
both the /+ and /_ eigenstates to form the TGBc p phase with either a fi or fn profile of 
the slab. 

Therefore, we conclude that the N*-TGBct transition is always split by the intermediate 
TGB 2q phase. At high b the period of oscillations of f±{x) becomes larger than its localiza- 
tion length, and the lowest eigenstate f+(x) has the single-peak profile of the TGBa phase. 
The N*-TGBA-TGBct-TGB 2q tetracritical point M corresponds therefore to the highest 
in b intersection of t± eigentemperatures. As follows from Fig. 3, the appearance of the 
TGB 2q phase becomes practically invisible when b decreases below the second intersection 
of t± at the point M\. We expect, therefore, that only the odd TGB 2q phase between Mo 
and Mi can be observed in reality. 

Equation ([0]) can be solved analytically in the two limit cases cii = and an ^> 1, when 
the operator Ti can be truncated to a more simple form. Although the case mi = does not 
correspond to the real situation of an > 1 we consider it first since it clarifies the qualitative 
structure of the phase diagram that conserves also at large an . 

B. Case cry = 

When a\\ =0, operator TC is a polynomial Schrodinger operator for the harmonic oscillator: 

U = a\ - 2a ± (-dl + b 2 x 2 ) + (-dl + b 2 x 2 ) 2 (14) 

Therefore, it has the same set of oscillator eigenfunctions: 

f n (x) = H n (Vbx)e~^ 2 (15) 

where H n are the Hermitian polynomials. The corresponding eigentemperatures are given 
by the equation: 

t n = 2a ± (2n + 1)6- (2n + I) 2 b 2 (16) 

The oscillations of t c2 (b) = maxt n (6) (Fig. 3a) are related to the quantum number n c2 of the 
lowest eigenstate that changes with b, unlike what occurs in the harmonic oscillator where 
the n = eigenlevel is always the lowest one. 
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When b > cr±/2, the n = eigenlevel does correspond to the upper critical temperature 
t c2 = 2a±b — b 2 . The block profile function has a Gaussian shape e~ x b ^ 2 in the TGBa phase. 
The slab width is calculated as in [|TJ: 

h = 2.2e l/2 /b 1 ' 2 (17) 

When b < a±/2, n c2 is given by <Tj_/26 — 1/2 rounded to the nearest integer. The oscillating 
upper critical temperature t C 2(b) tends to a\ when b vanishes as shown in Fig. 3a. The slab 
width is given by the width of the polynomial H nc2 (\/bx) as: 

h = crT/b (18) 
The block profile function f n ca(x) oscillates as cos(Ax + 7m c2 ) with period A = 2tt/ (2n c2 b) 1 ^ 2 

1 /2 

~ 2n/a/ . This corresponds to the slab of the TGB 2q phase which is the superposition of 
two SmC slabs with layers inclined at 6 = ±a]^ 2 . The parity of f n (x) alternates with n 
resulting in the different parity of f n< , 2 and f nc2 ±i eigenstates. So, the admixture of if nc2 ±i 
eigenf unctions with the TGB 2q profile function /„ c2 that occurs at some critical temperature 
below t c2 corresponds to the TGB 2q -TGBct transition. 

The tetracritical point M , where N*-TGBA-TGBct-TGBct phases meet is given by bo = 
o-j_/2, t = 3ai/4 

C. Case cr,| > 1 

In the limit <j\\ ^> 1 one can neglect b 2 x 2 in (|TT| ) and rewrite H as: 

n = (d 2 x + a ± f + a^x 4 (19) 



The eigentemperatures and eigenfunctions of Eq. ( |19D have in additional to (12) the scaling 
properties: 

t n (b,a h a ± ) = al-t c2 {ba\ / */a u l,l), (20) 
f n (x, 6, <t\\,<t±) = f n (aY 2 x, 6crJ /4 /a ± , 1, 1). 

So, t c2 (b, (T||, a±) at large cry is defined by the universal function t c2 (b, 1,1). We therefore 

use the scaled coordinates baV^/a±, tja\ for large cry to trace t c2 (b) that is obtained from 
numerical diagonalization of ([19]) (Fig. 3c). 

An important conclusion that followes from Fig. 3c is that the phase diagram for large cry 
posseses the same features as for cry = 0. It includes the domains of the TGBa and TGB 2q 
phases that precede the transition to the TGBct phase. The TGBa phase exists above the 
tetracritical point Mq and the TGB 2q phase is important in between the points Mq and M\. 
From Fig. 3c one locates the points M , Mi at ba^/a^ ~ 0.36 and ba}/ 4 /a± ~ 0.19. 

Consider now the transition to the TGBc phase at low field b 0.36o"x/o"^ 4 (i.e., below 

M ) and to the TGBa phase at high field b ^> 0.36<7_i_/cr^ 4 (i.e., above M ). Both the 
asymptotics for t c2 {b) are shown in Fig. 3c by the dashed lines. 
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Low fields: b <C 0.36cr_i_/cr^ 4 

At low field when the block width 7& diverges, one can assume that blocks have the structure 
of SmC slabs where the inclination of layers to the pitch axis is equal to its bulk value 
6*o = iarctancr^ 2 . Then the block's profile function f(x) has a high-frequency modulation 

, . l/2_ 

e ±la i_ x where + and — signs correspond to right and left tilting respectively. This factor 
can be excluded by shift: d x — » d x ± ia 1 / 2 in Eq.(|l9D. Next, we neglect d x in comparison 



with 2o~ 1 / 2 d x (justification of this approximation will be given below) and write 7i as: 

H= - (2(Ti /2 ^) 2 + C r||6 4 X 4 (21) 

Rescaling x as: ( = (6 4 cr||/4crj_) 1//6 x simplifies (|21"| ) to an unharmonic oscillator operator: 

n=(l6a 2 ± a ll ) 1/3 b^(-d 2 c + C i ) (22) 



Finally, t C 2(b) and the corresponding eigenfunction of ([H]) are given by: 



t c2 = a \ - 1.06 ■ (16^ cry) V3 6 4/3 (23) 



/ c2 (x) =e ± ^^ 1 ((6V||/4a ± ) 1 / 6 x 



where 1.06 and g\ (() are the lowest eigenvalue and eigenfunction of the unharmonic oscillator 
equation: (—9% + ( A )gi = 1.06g\. The block width l b is estimated as the characteristic width 
~ (4cr^/6 4 <T||) 1//6 of the function g x . More accurately, we can calculate 4 in the same way as 
was done in U for the TGBa phase, using the suitable Gaussian approximation for g\ (() 
as e -0,4 '' . Finally we get: 



J b = 3.2e 1 / 2 (o- x /&V||) 1 / 6 (24) 
At first glance, the eigenstate (^) is doubly degenerate with respect to a sign change in 

, . l/2_ 

the modulation phase e -l x . However, this is an artifact of the approximation: the term 

2 

d x we are neglecting lifts the degeneracy and splits the transition onto two. In particular, in 

— 2 

between the points M and M\ (Fig3c) d x favors the odd TGB 2q phase given by the linear 
combination (e 10 "^ x — e~ %a ^ x ) ■ g\ = 2igi sin o^x. The TGBct profile appears by means of 
an additional transition when the even combination ~ 2ig x cos cr ± x admixes. 

i 2 \ —2 

This effect, however, is as tiny (and negligible below Mi) as \d x ) ~ l/l h is smaller than 



2cr 1 / 2 <9 :r ^ ~ 2a 1 / 2 /lb, where lb is given by Eq.(|4]) and b 0.36crj_/cr^ 4 . The above estima- 
tion confirms the selfconsistency of the approximation we made. Condition b ^ 0.36a±/cr^ 4 

— 1/2 

is equivalent to the requirement that the wavelength of the transversal modulation A ~ a ± 
is smaller than the block width lb given by (p4]) . This clarifies the physical meaning of the 
approximation. 
High fields: b ^> 0.36crj_/(T^ 4 

Neglecting now a± in comparison with d 2 x we simplify Ti as: 
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n ~ sr x + a||6 4 x 4 (25) 

Then t c2 and the corresponding eigenfunction are given by: 

tc2 ~ _i. 39 . b 2 fc2 ft = g2 (al'h^x) (26) 

where 1.39 and g 2 (if) are the lowest eigenvalue and eigenfunction of the operator <9 4 + if. 
Unlike the low field case, the eigenstate fl25|) is nondegenerate and described by the even 
bell-shaped function g 2 , the lowest odd eigenstate being well separated from g 2 as shown in 
Fig. 3c. 

The eigenfunction g 2 fcrjj & 1//2 xJ gives the TGBa - like profile of the block. The block width 
lb is calculated using the results of U by noting that g 2 (tj) ~ e -°- 55r ? 2 : 

l b = 2.l8 1 l 2 /a\ / % 1 / 2 (27) 

/ — 2 \ —2 

The approximation we made is self-consistent since \d x ) ~ l/l b is indeed greater than cr_i_ 
when b ^> 0.36a±/a}/ 4 . Considering the truncated term h = 2a±d x + perturbation, 



one improves the asymptotic for t c2 (b). The correction to t c2 is given by (g 2 
27<jy //4 (j^6 + o\ where 7 = (g 2 d 2 g 2 \ — 0.5. Finally, we get: 
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t c2 ^al /4 a ± b-1.39al /2 b 2 (28) 
= -1.39 • (orj /4 b - 0.36a ± ) 2 + 0A8a 2 ± 

In fact, the asymptotic ( p8|) is valid also in the region of small negative <J±/b. It makes the 
analogous expressions in [114 , 1 5 1 more precise and explicitly determines all the coefficients. 



D. Stability of the TGB ct state. 

Having calculated the transition temperature and structure of the TGB ct phase we discuss 
now its stability, namely, whether this phase has the highest upper critical temperature t c2 
(or, alternatively, the highest upper critical field, b c2 ) among other TGB states and under 
which condition a transition from iV* to SmC* occurs via the intermediate TGBct phase 
rather than via the direct first order transition. 

To answer the first question, consider the competing TGBct and TGBc p phases that can be 
formed at the N*-SmC* transition. Dozov's estimation of the ratio of their upper critical 

fields 0], gives the factor 2 (C\\/2Dq$) = 2(2 

a \\) which is larger than one as long 



as CT[i > 1. Therefore, b c2 (TGBct) > b c2 (TGBc p ) and the TGBct phase is indeed more 
preferable. 

Our calculations are in agreement with Dozov's estimation. To show this, compare 
t c2 (TGB ct ) given by Eq.fl23|) with t c2 (TGB Cp ) calculated in |15j: 

t c2 (TGB Cp ) ~ u\ - 1.06 • (l + ay) 1/3 (l6ai<X||) 1/3 6 4/3 (29) 
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(We recalculated the numerical factor 1.06). The ratio of the upper critical fields for the 

TGBct and TGBc p phases (at given temperature) is expressed as (l + a\\\ ^ , which is larger 
than one. This is consistent with Dozov's estimation when an 3> 1, but gives a more precise 
value of the numerical factor. Although the starting GL functional (|]) is slightly different 
from that, used in |Tj| (see our remark after Eq.(fD) , this does not change the final result: 
one can show that in a relevant limit, a\\ ^> a±, expressions for b c2 {TGBct) calculated for 
both the functionals are the same. Note that standard perturbational analysis proves that 
the TGBct and TGB a phases are stable with respect to small parallel tilting of the layers 
when angle (p in Eq.(|9|) becomes nonzero. 

To answer the second question, note that transition via an intermediate TGB state is 
not preempted by the direct N*-Sm transition when the upper critical temperature t c2 for 
the TGB state is larger than the thermodynamical temperature t c of the direct transition. 
Recall first the calculations for the TGB a phase when a± < 0. Comparison of t c2 given by 
Eq.(|l~3"D (without (1 + 3<T||/4)6 2 ) with the N*-SmA transition temperature T n *a ©, which 
in dimensionless units is written as: 

t c = {gK 2 ) x l%/Dql (30) 
gives the following criterion for stability of the TGB a phase: 

k 2 > 1/V2 where k 2 = {gK 2 /2a\) l/2 /2Dql (31) 



When a± > the t c2 transition temperature for TGBct is given by Eq . (p5|) , and the 
N*-SmC* transition temperature T^*c* (§) in dimensionless units is written as: 

t c ~ a 2 ± - igI i 2)1 /\ l - a ± K 2 /2K 3 )b (32) 

The criterion t c2 > t c is always satisfied when b is small enough, since at small b, t c is 
linear in b and t c2 is proportional to 6 4,/3 . So, strictly speaking, one can always obtain the 
TGBct phase by preparation of the nearly racemic binary mixture of the left and right chiral 
molecules. The above consideration is valid only if t c (b) has a positive slope, that is, when 
cr± < 2K 3 /K 2 . This condition was used in [|l^,|nj as a criterion for stability of the TGBc 



state. It is not very likely, however, to find a system where this condition will not be the case 
since in conventional SmC* liquid crystals <jj_ < 1 and 2K 3 /K 2 ~ 4 — 6 [I5|. In contrast, the 



practical identification of this small-6 TGBct phase can be quite difficult because of the very 
large cholesteric pitch. We formulate the realistic criterion for the existence of the TGBct 
phase as a condition when it is stable in the vicinity of the critical point M . On the basis 
of the plot in Fig. 3c, one obtains that t c2 > t c near M if: 

k 2 > 0.5aJ /4 . (33) 



At realistic values of a\\ ~ 10 — 100 [24] this condition practically coincides with criterion for 
stability of the TGB a phase ([TIJ). The conditions ( |3"T| ) and (^) determine the position of 
the critical points Ma and Mc on the N*-SmA, C* transition line of Fig.l where the TGB a 
and TGBc phases first appear. 
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E. Resume 



We have revised the calculations of ]T5| for the N*-TGB^ct phase transition. The line of 
the upper critical temperature on phase diagram of Fig.l is reconstructed on the quantitative 
level and can be present in the relevant limit cry > 1 as follows: 

On the left side of the point M A that is defined by condition (|31|), the direct N*-SmA first 
order transition takes place. The critical temperature of the transition is given by flSOf) . 

Between the points Ma and M Q the transition to the TGB a phase occurs. The point M 
is placed at 0.3Qa± ~ a \\ A ^-> i ns ide the region a± > 0. The upper critical temperature t c2 of 
the transition is given by flT3| ) when —a±/b > 0.5 + 0.38<T||, and by (fZE|) when <r±/b varies 

from a small negative value to a±/b ~ a~y^ 4 /0.36. Both the asymptotics are matched in the 
intermediate region of negative cr±/b. 

In between the points M and M c the sequence of transitions N*-TGB 2q -TGB C t takes 
place, the location of the point M c being given by the condition (|33|). The TGB 2q phase 
exists in a small temperature interval below t c2 . It practically disappears on the right of 
the point Mi where a±/b > (X|j^ 4 /0.19. The upper critical temperature in between Mq and 
Mq is given by Eq. fl23|) . The junction of t c2 transition lines for TGB A and TGB 2q /TGB C t 
phases forms a kink in the tetracritical point M . 

On the right side of the point Mq the direct first order N*-SmC* transition occurs at the 
critical temperature t c , as is given by Eq.(|32|) 

We have also calculated the width of the TGB slab at T = T c2 with respect to different 
parameters b,a_i,a\\ (Eqns. (|T7|), flT8|), (^|) and fl2T|)). We summarize these results in Table 

I by the ratio h/ld = h/ld = ^b&/27r, resulting from the topological constraint ([[]). 



IV. ANALOGY WITH SPACE-MODULATED SUPERCONDUCTIVITY 

In this section we discuss a remarkable similarity between the TGB state and the mixed 
(vortex state) in superconductors of type II noted first in . We show that the TGBc phase 
corresponds to the mixed state in the "exotic" superconductors with space modulated order 
parameter. 

Note first that, due to de Gennes the iV — SmA transition is analogous to the 



superconducting phase transition because of similarity of the transformational properties of 
the order parameters: the space modulated function ip(r) = ip e qonr in SmA and the complex 
wave function \l/(r) in superconductor. The translation of ip(r) along the modulation vector 
qon (that is equivalent to multiplication of ip on a phase factor) corresponds to the gauge 
transformation in the superconductor. 

The long range twist of n(r) in iV* plays the role of the magnetic field destroying supercon- 
ductivity. The N*-SmA transition via the intermediate TGB a phase occurs with continuous 
Meissner-like expulsion of the twist of n. The rows of the screw dislocations resemble the 
Abrikosov vortex lattice. 

The SmC phase is characterized by the additional, transversal to n modulation of the 
order parameter resulted from the negative gradient terms in the CL model. Continuing 
the analogy with superconductivity one can tell that N-SmC transition corresponds to the 
transition to the superconducting state with nonuniform, space modulated order parameter. 
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The TGBc phase therefore should be an analog of the mixed (Abrikosov) state of modulated 
superconductor, providing that magnetic field is perpendicular to the modulation. 
On the quantitative level, the modulated superconductor in a magnetic field is described 
by the GL equation: 



tV» + </|*r* = C7j.(iV- Ay^ + D (zV- A)*<0 (34) 

where A is the vector potential of the field: B = curl A. The modulation of the supercon- 
ducting order parameter is provided by the gradient term with C± < 0. The upper critical 
field of this superconductor is calculated by neglecting the nonlinear term g\^\ \& and solv- 
ing the corresponding eigenproblem. It is easy to show that if A is chosen in the Landau 
gauge (0,Bx, 0), this procedure is exactly reduces to the case "mi = 0" we considered in 
Sect. IIIB. According to calculations of Sect. IIIB, this superconductor should have the 
oscillating upper critical field, shown in Fig. 3a. 
It is interesting to note that although the theory of space modulated superconductor was 
considered by Larkin and Ovchinnikov and Fulde and Ferrel |2T| more than thirty years 



ago, to our knowledge no direct experimental observation of this phase currently exists. 

Note finally another interesting analogy between the N*-TGBj\,ct phase transition and the 
B c2 transition in "unconventional" multicomponent superconductor UPt^ where the kink- 
like behavior of B C 2 was observed f22fl . The kink point in B c i in UPt^ is similar to the 
point Mo in the N*-TGBA,ct phase diagram: both points are provided by the intersection 
of eigenstates of the corresponding linearized GL equations. 
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TABLES 



TABLE I. Location of the tetracritical point Mq and the ratio h/ld in a different limit cases. 
Parameter e varies from when Kx^/Ki = to 1.5 when K\^/K2 is large. 





Location of Mq 






h/h a 




bo 


to 


TGB A , 


TGB 2q ,ct, 








b>b 


b<b 


G\\ <C 1 


0.5CT_L 


0.75a 2 ± 


0.8e 


a ± /2nb 


<7|| > 1 


0.36a ± /aJ /4 


0.43cri 




e(4a ± /ba ll ) 1 / 3 



a The ratio h/ld is given at N*-TGB transition line. At lower temperatures h/h increases rapidly. 
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FIG. 1. The phase diagram of TGB phases. The model parameters t, a± are controlled by the 
experimental conditions. We predict a new TGB<i q phase and penetration of the TGBa phase in 
the SmC* region where a± > 0. Dotted lines a, b and c present the N — SmA — SmC diagram in 
the nonchiral case. 
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FIG. 2. Structure of different TGB phases. 
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FIG. 3. Oscillating behavior of the lowest eigenstates of the lineariezed CL model that cor- 
respond the Cholesteric (N*) - TGBct transition as a function of the model parameters: t/a\, 
b/a± and cry. The concurrence between two nearly degenerate lowest eigenlevels (the minimal one 
corresponds to t C 2) leads to spliting of the transition by appearance of either the TGB\ phase or 
the intermediate TGE>2q phase. Dashed lines in (c) present the asymptotics calculated in the text. 
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